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Abstract 

Based on the shearlet transform we present a general construction of 
continuous tight frames for L 2 (R 2 ) from any sufficiently smooth function 
with anisotropic moments. This includes for example compactly sup- 
ported systems, piecewise polynomial systems, or both. From our earlier 
results in _5 it follows that these systems enjoy the same desirable ap- 
proximation properties for directional data as the previous bandlimited 
and very specific constructions due to Kutyniok and Labate; [8]. We also 
show that the representation formulas we derive are in a sense optimal for 
the shearlet transform. 
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1 Introduction 

The purpose of this short paper is to give a construction of a system of bivariate 
functions which has the following desirable properties: 

directionality. The geometry of the set of singularities of a tempered 
distribution / can be accurately described in terms of the interaction be- 
tween / and the elements of the system. 

tightness. The system forms a tight frame of L 2 (R 2 ). 

locality. The representation is local. By this we mean that the represen- 
tation can also be interpreted as a representation with respect to a non 
tight frame and its dual frame such that both of these frames only consist 
of compactly supported functions. 
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The construction is based on the shearlet transform which has been introduced 
in [S] and has become popular in Computational Harmonic Analysis for its 
ability to sparsely represent bivariate functions. A related important result 
is that the coefficients of a general tempered distribution with respect to this 
transform exactly characterize the Wavefront Set (see e.g. [7] for the definition) 
of this distribution [5]. 

Notation. We shall use the symbol | ■ | indiscriminately for the absolute value 
on R, R 2 ,C and C 2 . We usually denote vectors in R 2 by x,t,£,u> and their 

elements by Xx,X2,t±,t2, In general it should always be clear to which 

space a variable belongs. The symbol || • || is reserved for various function space 
and operator norms. For two vectors s,f€l 2 we denote by st their Euclidean 
inner product. We use the symbol / < g for two functions /, g if there exists 
a constant C such that f(x) < Cg(x) for large values of x. We will often 
speak of frames. By this we mean continuous frames as defined in pQ. We 
define f(uj) — J f(x) exjp(27riujx)dx to be the Fourier transform for a function 
/ G L 1 n L 2 and continuously extend this notion to tempered distributions. 

1.1 Shear lets 

We start by defining what a shearlet is and what the shearlet transform is: 

Definition 1.1. A function ip g L 2 (R 2 ) is called shearlet if it possesses M > 1 
vanishing moments in x\-direction, meaning that 



'I, doj < oo. 

2 |wi| 



2M 



Let f G L 2 (R 2 ). The shearlet transform of f with respect to a shearlet ip maps 
f to 

SH^f(a,s,t) := (f,ip ast ), 



whe 



tp ast (x) := ^( , —^r, 



a 

The following reproduction formula holds [3] : 
Theorem 1.2. For all f G i 2 (R 2 ) and tp a shearlet 



f(x) = / / / S / H % i,{a,s,t)%l) ast (x)a dadsdt, (1) 
Jm 2 Jr jr + 

where equality is understood in a weak sense. 

The shearlet transform captures local, scale- and directional information via 
the parameters t, a, s respectively. A significant drawback of this representation 
is the fact that the directional parameter runs over the non compact set R. Also 
it is easy to see that the distribution of directions becomes infinitely dense as 
s grows. These problems led to the construction of shear lets on the cone [5]. 
The idea is to restrict the shear parameter s to a compact interval. Since this 
only allows to caption a certain subset of all possible directions, the function 
/ is split into / = Pf + P v f, where P is the frequency projection onto the 
cone with slope < 1 and only Pf is analyzed with the shearlet ip while P v f 
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is analyzed with i\) v (x\,X2) '■— VK^j^i)- For very specific choices of ip Labate 
and Kutyniok proved a representation formula 



\\f\\l = / \(f,T t W)\ 2 dt+ [ff \Sn^Pf(a,s,t)\ 2 a- 3 dadsdt 

+ [ [ [ \SU r P v f{a,s,t)\ 2 a- 3 dadsdt, (2) 
Jr 2 J -2 

T t denoting the translation operator by t. The main weakness of this decompo- 
sition is its lack of locality: indeed, first of all the need to perform the frequency 
projection P to f destroys any locality. But also the functions ip which have 
been considered in [8] are very specific bandlimited functions which do not 
have compact support. As a matter of fact no useful local representation via 
compactly supported functions which is able to capture directional smoothness 
properties has been found to date. The present paper aims at providing a step 
towards finding such a representation using two crucial observations: First, in 
[5] we were able to show that the description of directional smoothness via the 
decay rate of the shearlet coefficients for a — > essentially works for any func- 
tion which is sufficiently smooth and has sufficiently many vanishing moments 
in the first direction, hence also for compactly supported functions. The second 
observation is that actually a full frequency projection P is not necessary to ar- 
rive at a useful representation similar to ©. Instead of the operator P we will 
use a 'localized frequency projection' which is given by Fourier multiplication 
with a function pq to be defined later. 

Our main result Theorem 12.41 will prove a representation formula similar to 
@ where ip is allowed to be compactly supported and the frequency projection 
is replaced with a local variant. We also show that this local variant of the 
frequency projection is the best we can do - without it, no useful continuous 
tight frames (or continuous frames with a structured dual) can be constructed 
within the scope of the shearlet transform. This is shown in Theorem [ 

2 The construction 

The goal of this section is to derive a representation formula 

(■1 r \ 



/111 = Tril I I \(f,qo*ipa S t)\ 2 a- 3 dadsdt 
W Jr 2 j -i Jo 

\{f, qi *^ ast )\ 2 a- 3 dadsdt 
\(f,T t v)\ 2 dt) (3) 
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for L 2 (]R 2 )-functions / and with some localized frequency projections qo,qi and 
a window function ip to be defined later. In order to guarantee that the shear let- 
part contains the high-frequency part of / and the rest contains low frequencies, 
it is necessary to ensure that the window function ip in this formula is sufficiently 
smooth. 
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2.1 Ingredients 

Here we first state the definitions and assumptions that we use in the construc- 
tion. Then we collect some auxiliary results which we later combine to prove our 
main results Theorems 12.41 and 12.61 The large part of the results will concern 
the construction of a useful window function and to ensure its smoothness. We 
say that a bivariate function / has Fourier decay of order Li in the i-th variable 
(i€{l,2})if 

We start with a shearlet tp which has M vanishing moments in si-direction 
and Fourier decay of order L\ in the first variable. It is clear from the definition 
of vanishing moments that t/j — (^§-) with 8 g L 2 (R 2 ). We assume that 
9 has Fourier decay of order L 2 in the second variable so that the following 
relation holds: 



We also set 



2M- 1/2 > L 2 > M > 1. (4) 



N :=2min(L 2 - M,Li), (5) 
C* ■= i *f (6) 
A^(0~ I I ma^a^ib-s^fa-^dads. (7) 



and 



-2 Jo 

We define functions (fo, ipi via 

l^o(0| 2 = C 4 , - A^(0 and |&(0| a = C«, - (£), (8) 

where 

We write \C for the indicator function of the cone C = {£ = (616) G K 2 : 
161 16 I}- Finally, we pick a smooth and compactly supported bump function 
<I> with $(0) = 1 and define functions po, p\ via 

Pa = $ * xc and pi = 1 - p . (9) 
Clearly, po and p\ are both compactly supported tempered distributions. 
Lemma 2.1. We have 

\Po(0\ < \t\- N M 161/161 > 3/2 

IM0I £ ier N M 161/161 > 3/2 (10) 

Proof. Assume that £ = te with e a unit vector with | ei | / 1 62 1 > 3/2 and < > 0. 
There exists a uniform S > such that for all r] with |r?| < <5£ we have £ — -q G C c , 
and hence xc(£ — ^) = 0. It follows that we can write 

|Po(OI = l/ Xc(Z - v)Mv)dri\ < f \Hv)\dr) 

JR 2 J\r)\>St 
< t -N = ^-N 
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if $ is sufficiently smooth. On the other hand, let £ = te with e a unit vector 
with 1 62 1 / 1 ei | > 3/2 and t > 0. There exists a uniform 6 > such that for all 77 
with \rj\ < St we have £ — 77 6 C and hence xc(£ — ^) = !• Now we can estimate 



|ft (01 = |l-0o(0l = ly 2 Hr))(l-Xc(Z-v))dr)\ 

Hv)(i-xc(t-v))dv\<t- N = \a\- N 



\ v \>6t 

again for $ smooth. Note that in the first equality we have used that $(0) = f . 
This proves the statement. □ 

Lemma 2.2. We have 

\M0\ 2 < \t\- N for < 3/2 

l^i(0| 2 < ft-" for 161/161 < 3/2 (ff) 

Proof. This follows from [5J Lemma 4.7]. For the convenience of the reader we 
present a proof here as well. We only prove the assertion for ipo since the proof 
of the corresponding statement for ipi is the same. By definition we have 

|^o(6l 2 = (/ |^(a6,^(6-s6))| 2 a- 3/2 da<fs 

JaGR, |s|>2 

I^(a6,^(6-s6))| 2 a~ 3/2 ^). 

a>l, |s|<2 

We start by estimating the second integral using the Fourier decay in the first 
variable: 

f |^(a6,v^(6 -s^))\ 2 a- 3 / 2 dads < 4 f (a|6l) _2Ll a^da 

Ja>l, \s\<2 Ja>l 

< i6r 2ii <ier 2Li 

for all £ in the cone with slope 3/2. To estimate the other term we need the 
moment condition and the decay in the second variable. We write -0(6 = 
and £ = (£i,r£i), \r\ < 3/2. We start with the high frequency part: 

|^(a6,^(6-s6))| 2 a" 3/2 da 

a<l, |s|>2 



< 



/ |a6| 2M |^(a6,V^(6-s6))| 2 a" 3/2 ^ 

Ja<l, \s\>2 

I |a6| 2M |v^(6-.s6)r 2L2 «- 3/2 ^ 

Ja<l, \s\>2 

f a 2M-L 2 -3/2| 6 |2M-2L 2 | r _ s |-2L 2dads 
Ja<l. \s\>2 



'a<l, |s|>2 

we now use that \r — s\ is always strictly away from zero. By assumption 
L2 = 2M — 1/2 — e for some e > 0. Hence we can estimate further 

_ = | a |-2(i 2 -M) f a - 1+e \r-s\- 2L »dads<\^ L '-^. 

Ja<l, \s\>2 
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The low-frequency part can simply be estimated as follows: 

|$(a&,Vo(6 - s^))\ 2 a- 3/2 dads 



I 

J a 



a>l, |s|>2 

< |Ci|- 2L2 / a- 3 / 2 - L *\r-s\- 2L *dad S 

Ja>l, \s\>2 

< ir 2i2 - 

Putting these estimates together proves the statement. □ 

Lemma 2.3. Assume that %p is compactly supported with support in the ball 
Ba '■= {£, ■ |£| < A}. Then the tempered distributions (|<£>j| ) v , i = 0, 1 are both 
of compact support with support in the ball 2 y3 + \/%Ba ■ 

Proof. Since the inverse Fourier transform of the constant function is a 
Dirac, this follows if we can establish that the functions and A^* are Fourier 
transforms of distributions of compact support. We show this only for A^, the 
other case being similar. In what follows we will use the notation f~(x) := 
f(—x) for a function /. Since ip has M anisotropic moments we can write ip = 
(~5x^) M ^ or some ^ € L 2 (M. 2 ) with the same support as ip. Consider the function 
©(£) := l^(£)| 2 - It is easy to see that this is the Fourier transform of the so-called 
Autocorrellation function 6*8~ of 9 which is compactly supported with support 
in Bia- Therefore, by the theorem of Paley- Wiener, the function O possesses 
an analytic extension to C 2 which we will henceforth call F. Furthermore, by 
the same theorem F is of exponential type: 

\F(0\ <(l + |CI) K exp(2A|3CI), (12) 

where ( e C 2 and K € N. We now consider the analytic extension T of the 
function A^, which is given by 

r(C) = f / Ki) 2 ^Ki»« 1/2 (C2 - s^a^dads, c G c 2 . 



Since M > 1 by (jU the above integral is locally integrable which implies that 
r is actually an entire function. It is also of exponential type: Writing M as := 

a short computation reveals that 



-so 1 ' 2 a 1 ' 2 



\M as \\ LHC ^ L2(c2) < a^(i + s l + ( S 2 + S l)V*f 2 = : a ^C{s). 



Now we estimate 

p2 



A^(OI = 1/ / Ki) 2M ^(A^ s C)^ 3/2 ^ s |<|C| 2M (l + |M as C|) if exp(|M as 3CI) 
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< sup (1 + \(\) 2M+K exp(2AC(s)\-SC\) < (1 + \C\) 2M+K exp(2AC(2)|3C|). 

s£[-2,2] 

By the Theorem of Paley- Wiener-Schwartz 7 it follows that A^, is of compact 
support with support in ^ 2 y/ 3+y/E A' ^ 
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2.2 Main Result 

We are now ready to prove our main result, the local representation formula in 
Theorem 12 .41 which is similar to (0) only with local frequency projections (given 
by convolution with Po,Pi) and with possibly compactly supported shearlets. 
In addition, in Theorem 12.61 we show that in a way this is the simplest repre- 
sentation that one can achieve with shearlets - the local frequency projections 
are necessary in order to wind up with useful systems. 

From now on we shall assume that < $(£) < 1 for all (el 2 and define 

Furthermore, we define 

m ■■= (M0lw(0l a +M0M0l 2 ) 1/2 - 

Observe that by the positivity assumption above the radicands in the previous 
definitions are nonnegative and real. 

Theorem 2.4. We have the representation formulas 

ll/lla = 7r{ I I I' \(f,qo*i>ast)\ 2 a- 3 dadsdt 

L'V Jm. 2 J -2 JO 
2 ,1 

/ \(f, qi *r ast )\ 2 a- 3 dadsdt 

-2 JO 

\(f,T t v)\ 2 dt) (13) 



an 



d 



f = f hi9h + f l0W := ±- ( / / [ (f, ^ast) P o * iPasta-'dadsdt 

(■2 r l 

(f^ast)Pi * ^asta~ 3 dadsdt) 
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i(/ (f,T m )po*T t v dt+ I (f,T t <pi)pi*T t <pidt). (14) 



The function fl° w satisfies 



(f low r(o < \z\- N (is) 



for any f. We have the following locality property: If i\) is compactly supported 
in Ba 
t + B 



in Ba and $ has support in Bb, then f^ ow (t) only depends on f restricted to 



2^/3+V5A+B' 

Proof. We first prove (fT4")l . Taking the Fourier transform of the right hand side 
yields 

i-(p (£)(A^) + \M0\ 2 ) +pi(0(Av."(0 + |ft(Ol a ))/(0 = 
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The proof of (fT3"|) is similar. We prove the equation (fT5j) : Again taking the 
Fourier transform of f^ ow gives 

^(MOlMOl 2 +pi(0\MO\ 2 )m)- 

Now, the desired estimate follows from Lemmas 12.11 and 12.21 The last statement 
follows from the observation that f^ ow can be written as 

/ low =^-/*(po*(|^o| 2 ) v +W*(l<?i| 2 ) v ) 



together with Lemma \2. 31 □ 

Remark 2.5. While we could show that the functions tpj * ifi and pi, i = 1,2 
are compactly supported, it would be desirable to show that also the functions 
(pi, qi, i = 1,2 have compact support. We currently do not know how to do this. 
By a result of Boas and Kac Jlf, in the univariate case there always exists for 
any compactly supported function g with positive Fourier transform a compactly 
supported function f with f * / _ = g. However, in dimensions > 2 this holds 
no longer true and things become considerably more difficult. 

The previous theorem for the first time gives a completely local represen- 
tation for square integrable functions which also allows to handle directional 
phenomena efficiently: by the results of [5] it follows that the decay rate of the 
coefficients (f,i/j a st) for a ^ accurately describes the microlocal smoothness 
of / at t in the direction with slope s. 

It is interesting to ask if the frequency projections given by convolution with 
Po,Pi are really necessary, or in other words if it is possible to construct tight 
frame systems (T t (p) teM 2 U {ipast)ae[o,i], s e[-2,2],tm 2 u (C S i)a£[o,i], s e[-2,2],teK 2 for 
L 2 (R 2 ). We show that this is actually impossible, meaning that in a sense 
Theorem 12.41 is the best we can do. 

We remark that for any shear let if), the constant can also be computed 

as 

/oo />oo 
/ |?A«i,a 1/2 (6 - s^))\ 2 a~ 3 / 2 dads, 
-oo JO 

as a short computation reveals. This fact is related to the inherent group struc- 
ture of the shear let transform [4]. 

Theorem 2.6. Assume that tp — (j>§^) M $ is a shearlet with M > 1 vanishing 

moments in x\-direction such that either M > 1 or M = 1 and 9 G i°°(]R 2 ). 
Furthermore we assume and L\ > 0, Li > M with L\,Li defined as in Lemma 
\2."A Then there does not exist a window function (p such that 

lim = 

and such that the system (T t ^) teR 2U(V'a S t)Qe[o,i],se[-2,2],teR 2 U(V'a S t)ae[o,i], s e[-2,2],iGi 
constitutes a tight frame for L 2 (R 2 ) , which means that a representation formula 



1 



l/H 2 = ^( / / / \{f^ast)\ 2 a- 3 dadsdt 



2 Jo 

1 \(f,r a st)\ 2 a- 3 dadsdt+ f \(f,T tV )\ 2 dt) (16) 
Jr 2 
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holds for all f £ L 2 (M. 2 ) and some constant C . 

Proof. In terms of the Fourier transform, (|16[) translates to 

Mo + A^(o + |0(Ol a = c. 

If we assume that ip has Fourier decay lim^i^^ <p(£) — this would imply that 

lim (A^(0+A^.(0) =C. 

|£|-s-oo 

Lemma [2?2l implies that lim^i^oo A,/,(£) = lim|£|_ yoo A.0^(£) = C^, for | < J|4 < 

|. It follows that C = 2C^. On the other hand, using the moment condition 

$ = £f§ for M = 1 and 6 G L°°(]R 2 ), we have the following estimate for A^, 
and £ in the strip Ss := {£ : |£i| < 5}: 

IMOI < ^ 2 / / a^lflllooo-^dodfl < SII^IU^ 2 . 

J -2 JO 

If we assume that M > 1, we can write ^>(£) = £i* _: V(£) where fj. is still a 
shear let. We get a similar estimate as above: 

J-2 JO 

At any rate, by choosing S appropriately small, this implies that for t; £ Sg we 
have 

C=lim (A*(0 + A^(0) <C, 

K|->oo 

which gives a contradiction. □ 

It is easy to extend this argument to show that there do not exist shearlet 
frames such that there exists a dual frame which also has the structure of a 
shearlet system, see [1] for more information on frames in general and [6] for 
more information on shearlet frames in particular. 

3 Concluding Remarks 

In future work we would like to address the problem of constructing continuous 
tight frames for the so-called 'Hart Smith Transform' [3J [IT] where the shear 
operation is replaced with a rotation. We think that in this case the results 
might become simpler. One reason for this is that in this case no (smoothed) 
projection onto a frequency cone is needed. In view of constructing discrete 
tight frames we think that a simple discretization of a continuous tight frame 
will not work for non-bandlimited shearlets. The approach that we are currently 
pursuing in this direction is to construct so-called Shearlet MRA's via specific 
scaling functions and to try to generalize the 'unitary extension principle' of 
Ron and Shen [TU] to the shearlet setting [B]. 
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